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Abstract 

It was shown that the Runge-Lenz vector for a hydrogen atom is equivalent to the 
raising and lowering operators derived from the factorization of radial Schrodinger 
equation. Similar situation exists for an isotropic harmonic oscillator. It seems that 
there may exist intimate relation between the closeness of classical orbits and the 
factorization of radial Schrodinger equation. Some discussion was made about the 
factorization of a ID Schrodinger equation. 



PACS number(s), 03.65.-w, 03.65.Fd, 03.65.Ge 



I. INTRODUCTION 



The factorization method of Schrodinger '^'^1 was extended to address the radial 
Schrodinger equation of a particle in a central potential V{ry^~^^. It was shown that 
only for two kinds of potential, the Coulomb potential and isotropic harmonic oscillator, 
the radial Schrodinger equation can be factorized, and directly from the factorization of 
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radial equation, two kinds of raising and lowering operators, A±{1) and B±{1), can be 
derived for an isotropic harmonic oscillator, but only one, A±{1), for a hydrogen atom. 
This reminds us of the famous Bertrand theorem'^l in classical mechanics, which says: 
"the only central forces that result in close orbits for all bound particles are the inverse 
square law and Hooke's law" . It is guessed that the factorizabihty of radial Schrodinger 
equation may be intimately connected with the closeness of classical orbits. 

In classical mechanics, the maximum number of functional independent conserved 
quantities of a closed Hamiltonian system with N degrees of freedom is 2A'" — l^^l For a 
system with independent conserved quantities no less than N is called integrable^^l . An 
integrable classical system with A^+A independent conserved quantities (0<A<A"— l)is 
called A-fold degenerate, and there exist A hnear relations with integer coefficients between 
the A^ frequencies u!i{i = 1,2, N) oi the systeml^°l . A classical system for A = A^ — 1 is 
called a completely degenerate system, and there remains only one independent frequency. 
For example, it is well-known that for a particle in a central potential V{r), apart from 
the Hamiltonian, the angular momentum 1 is also conserved, and the particle in a general 
central potential V{r) is 1-fold degenerate and moves in a plane, but the planar orbits are 
in general not closed. For a particle in an attractive Coulomb potential {V{r) = —n/r), it 
was demonstrated [^^1 that there exists an additional conserved quantity, the Runge-Lenz 
vector, (m — K — l),a = p x 1 — r/r, which guarantees the closeness of a Kepler orbit. In 
fact, the direction of a is just that of the major axis of the elliptic orbit and the magnitude 
of a is the eccentricity. The existence of Runge-Lenz vector implies that the Coulomb 
potential has a higher dynamical symmetry SO 4, than its geometric symmetry SO^ 
However, it is seen that a • I = 0, and a? — 2H]? -\- 1, so the number of independent 
conserved quantities is 5, and the hydrogen atom is a completely degenerate system. 
Similar situation exists for an isotropic harmonic oscillator. 
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In Sect. II it will be shown that for a 2-dimensional(2D) hydrogen atom, the Runge- 
Lenz vector is equivalent to the raising and lowering operators derived from the factoriza- 
tion of radial Schrodinger equation. In Sect. Ill, it will be shown that, for a 3D hydrogen 
atom, from the Runge-Lenz vector a and angular momentum I, one can construct three 
kinds of raising and lowering operators {An — 0, AZ = ±l,Am = 0, ±1), which are 
equivalent to the raising and lowering operators derived from the factorization of radial 
Schrodinger equation. In Sect. IV we will address isotropic harmonic oscillators. In terms 
of the two kinds of raising and lowering operators one can construct the conserved quan- 
tities characterizing elliptic orbits. Finally, in Sect. V there is some discussion about the 
factorizability of the Schrodinger equation for a ID system, which has been investigated 
extensively in super symmetric quantum mechanics. 

II. 2D HYDROGEN ATOM 

(A) Runge-Lenz vector and a simple algebraic approach to the 
eigenvalue problem 



For a 2D hydrogen atom, the quantum version of the Runge-Lenz vector reads 
{h — m — e — 1) 

a = ^(p X I - I X p) - Bp = p X I - ip - Cp, (1) 
where p = Pxi + Pyj, 1 = i^k = {xpy — ypx)^, P — y/x^~+y^- It is easily shown that 

[h: cix] = iay, [1^, ay] = -iax, [a^, ay] = -i2Hl^. (2) 

In the subspace spanned by the bound states with eigenenergy E[E < 0),ax and ay may 
be replaced by — —l/2Eaxi — y^— l/2£'ay. Let A^ — Iz, then 

[Aa,Ap] ^ iSafS-yA^, (3) 
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i.e. Ax, Ay and A^ constitute the SO3 Lie algebra, and the eigenvalue of = — | — 2^ is 
mo{mo + 1), mo = 0, 1, 2, • • •. Therefore, the energy eigenvalue of a 2D hydrogen atom is 

E^En^ -l/2n\ n = (mo + 1/2) = 1/2, 3/2, 5/2, • • • . (4) 

Alternatively, defining a± = a^; ± ia^, it is easily shown that 

[Z^,a±]=±a±, (5) 

a_a+ = y(2Z, + 1)2 + 1, (6) 

i.e. a± are just the raising and lowering operators of the magnetic quantum number m. 
Assume \Em) is the eigenstate of {H,lz) with eigenvalues {E,m), then a±\Em) are also 
the eigenstates of H with energy E and the eigenstates of Iz with eigenvalues (m ± 1). 
For a given energy eigenvalue the allowed |m| must have an upper limit mo, and 
a+|£^mo) = 0. Hence a_a+|£'mo) = 0. Using (6), we get f (2mo + 1)^ = 1, which 
is just (4). The degenerate states belonging to En may be expressed as a^|£^„mo). A; ~ 
0, 1, • • • , 2mo, and the degeneracy is 2n — (2mo + 1) = 1, 3, 5, • • • . 

(B) Raising and lowering operators derived from factorization 



The energy eigenstate of a 2D hydrogen atom may be chosen as the simultaneous 
eigenstate of {H, 1^), i.e. ^(p, 0) — ^'^'^Xm{p)l \fp-, m — 0, ±1, ±2, • • •, and Xm{p) satisfies 
H(m)xm{p) = Exm{p), or 

D{m)xm{p) = A^Xm(p), >^m = -2-B, (7) 
D{m) = -2H{m) = — ^ + -. 



Directly from the factorization one may derive the raising and lowering operators'^', 

, , , d m + 1/2 1 
A+{m) = '— + 



dp p m + 1/2' 

, , , d m — 1/2 1 , , 

A_{m) = - + ^ — , 8 

dp p m — 1/2 

whose selection rules are /S.E — 0(An = 0) and Am = ±1. Using A±{'m), the energy 
eigenvalues and eigenstates can also be easily obtained. From (7) it is seen that E depends 
only on the absolute value of m. Using Hellmann-Feynman theorem 

dE _mm)_^Jn^^^^ 



d\m\ d\m\ p^ 
i.e. E increases monotonically with |m|. Thus, m = for the ground state. On the other 
hand, for a given E, \m\ must have an upper limit, say, mo(mo > 0) . Then 

^.(".o)x™«W = (|; - + = 0- (10) 

So Xmo(p) ~ p'"o+i/2e-p/(mo+i/2)_ Yyoui (10) we have A_(mo + l)A+(mo)Xmo(p) = 0, and 
using (7) and (8), we get 

[L>(mo) - l/(mo + ll2f]xmM = [-'^E^, - + ll2f]xmM = 0. (11) 

Thus we get E^^ — — l/2(mo + l/2)^. Using >l_(mo), A_(mo — 1), • • •, successively operat- 
ing on Xmo(p)) may get all the degenerate eigenstates belonging to Em^, Xm(),m(p), ?7t. = 
mo, mo— 1, • • • , —mo. Alternatively, the degenerate eigenstates belonging to En = — l/2n^, (n 
1/2, 3/2, 5/2, • • •) may be denoted by Xn,m{p), \'^\ — n—l/2,n — 3/2, •••,1,0, and it can 
be shown that Xn,m{p) ~ pl'"l"'"^/^e~''/"'F(— n^, 2|m| + 1,^), where F is the confluent 
hypergeometric function, and = (n — 1/2) — |m| = 0, 1, 2, • • • , (n — 1/2). 

(C) Equivalence of the Runge-Lenz vector and the raising and 
lowering operators 
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In polar coordinates the Runge-Lenz vector (1) can be expressed as a± = ± lay 

Operating on the eigenfunction e^'^'^Rmip), the raising and lowering of m in the angular 
function are accomplished by e^*'^, and a± are equivalent to operators a±{m) operating 
on Rm{p), 

«.(m) = [(Tm^ + _)--- ± -- 1] 

,1 , d m(m± 1/2) , , 

Apart from a trivial constant factor, a±{m) may be expressed as 

dm 1 

a+{m) = h 



dp p m+1/2' 

, . d m 1 ,^ 

a4m) = — + — , 14 

op p m — 1/2 

The operators a±{m) operating on Rm{p) may be replaced by A±[m) operating on 
Xm{p) = y/pRm{p) 

^ , , d m + 1/2 , 1 
A+{m) = — + 



dp p m+1/2 

. / N <9 m-1/2 1 

A(m) = 7^ + -. (15 

^ ' dp p m-1/2 ^ ^ 

which are just the raising and lowering operators (8) derived directly from the factorization 
of the Schrodinger equation. 



III. 3D HYDROGEN ATOM 



Now we address the Runge-Lenz vector for a 3D hydrogen atom, 

a = pxl — ip — r/r. 
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(16) 



Defining a± = ± iay, we get 

In the spherical coordinate system a± may be expressed as 

d 

a± = ±— [cos^/± -sin^e^*'^(/^ ± 1)] 

:F^[sin^^Z± ± cosei±{h ± 1) T sin^e^^'^Z;,(Z^ ± 1)] 
-sin^e^^'^. 



Similarly, 



= ^[^sin^(e*^/_ - e-^'^/+) - cos^] 

+-{sin^— + -[cos^(Z+Z_ + - sin^(e''^Z^Z_ + e-'%Z+)]} 
— cos 9. 



(18) and (19) operating on the simultaneous eigenfunction of (if, 1^, Z^), ^'„;^(r, 6*, 
Rni(r)Yim(9,(l)), we get 

a±^nlm — ±^-RnK'^)[(^ + l)^/,±m^+l,m±l + ^C^i-l,-(±m+l)^-l,m±l] 

^^Rnl{r)[l{l + l)(iz,±m^i+l,m±l " + l)di-i-(±rn+l)yi-l,m±l] 
±Rnl{r)[di^±mYi+i^rn±l — dl-l,-{±m+l)Yl-l,m±l\, 
O'z'^nlm = —-^Rnl{'f')[{l+^)ci,mYl+l,m — lCl-l,mYl-l,m\ 

+ ^Rnl{r)[l{l + l)Q,„yj+i,^ + l{l + l)ci-i^rnyi-\,m\ 
— Rnl{f)[ci,mYl+l,m + Q-l,m^i-l,m] , 



where 



Using (20) and (21), we may get 



ar r I + 1 



nlm 



r d 

+t + 

, d 

|. d 

,d 
+^d-r + 



— - j]Rniir)lil - m - 1)q_i,^1^_i,„, 

]Rni{r){l + l){l-m+ l)di^myi+i,m+i 



l + l 

-j]Rni{r)l{l + m)(ii_i _(^+i)li_i,„,+i, 

+ J^]Rnl{r){l + 1)(Z + m + l)di,^mYl+l,m-l 

]Rni{T)l{l — rn)di-i,m-iYi — l,m— 1- 



(23) 



(24) 



(25) 



Now, we may define the operators S± constructed from tlie Runge-Lenz vector and angular 
momentum operator. 



5"+ = Z_a+ + (^2 - i + 1)0^, 



(26) 



It is easily verified that 



S+-^nlm = {-^-- + -^)Rnl{r){l + l){2l + l)ci,mYl+i,m{e,(l>) 

ar r / + 1 

oc a+(/)i?„;(r)F/+i,m(6',(/)), 

= {^ + ^-^-\)Rni{r)l{2l + l)ci_^,rnYi-i,^{eA) 
dr T I 

oc a-{l)Rni{r)Yi-i,m{0-,(t))-, 



(27) 



where 



, d I 1 . , (i Z + 1 1. 



(28) 



which are equivalent to the angular momentum raising and lowering operators A±{1) 
derived from the factorization of radial Schrodinger equation for a 3D hydrogen atom^^l , 



^ ,,, . d l + l 1 X . /,v , d I 1, 



(29) 
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which operate on the radial wavefunction x„;(r) = rRni{r). It is seen that the effect 
of is to increase (decrease) the angular momentum Z by 1, but keep the energy 

and magnetic quantum number m unchanged. Therefore, to clearly indicate the selection 
rules, S+ and S- may be relabelled as 

S-^S{n,li,m) = /_a+ + (/, + / + 2)a,. (30) 

Similarly, using the Runge-Lenz vector and angular momentum operator, one may con- 
struct the other two kinds of raising and lowering operators 





= {Iz + l - l)a+ - l+O'z, 




S{n,l i,m t) 


= {h-l - 2)a+ - l+ttz, 


(31) 


S{n,ltmi) 


= {lz — l + l)a- - l-ciz, 




S{n,l i,mi) 


= {lz + l + 2)a_ - La^. 


(32) 



IV. ISOTROPIC HARMONIC OSCILLATORS 



It is well-known that an nD isotropic harmonic oscillator has the dynamical sym- 
metry SUn- For a 3D isotropic harmonic oscillator, apart from the Hamiltonian H and 
angular momentum I, there exist five additional conserved quantities which constitute a 
quadruple tensor 

Qxy — ■^y ~l~ PxPyi Qyz ~l~ PyPzi Qzx — ~l~ PzPxi 

Qi = l[i^'-y') + ipl-pl)], (33) 

Qo = ^[{x' + y' - 2z') + {pI + pI - 2pl)]. 
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It can be shown that there exist four relations among the nine conserved quantities, so 
a 3D isotropic harmonic oscillator is also a completely degenerate system and moves, in 
general, along an elliptic orbit and the direction of the semi-axes and eccentricity are 
characterized by the quadruple tensor. It has been shown'^J that the radial Schrodinger 
equation of an isotropic harmonic oscillator as well as a hydrogen atom can be factorized. 
Nevertheless, it was noted that for an isotropic harmonic oscillator, two (rather than 
one) kinds of raising and lowering operators can be derived from factorization, and these 
operators themselves are not conserved quantities (unlike the situation of the hydrogen 
atom). However, it can be shown that in terms of the two kinds of raising and lowering 
operators one can construct the conserved quantities characterizing an elliptic orbit. 

For simplicity, we take a 2D isotropic harmonic oscillator as example. A 2D isotropic 
harmonic oscillator has the dynamical symmetry SU2, which is locally isomorphic to SO^ 
(the dynamical symmetry of a 2D hydrogen atom). For a 2D isotropic harmonic oscillator 
in classical mechanics, one may construct two conserved quantities, Qxy = xy + PxVy ^-iid 
Qi = — y^) + {pi — Py)]-! intersection angle 7 of the major axis with the x 

axis is determined by tan 27 = Qxy/Qi-, ^-nd the eccentricityoc [Q^y + <5i]^- 

In ref. [6], it was shown that from the factorization of radial Schrodinger equation 
for a 2D isotropic oscillator, two kinds of raising and lowering operators, A± and B±, 
operating on the radial wave function Xm(p), can be derived 




(34) 



which may be replaced by a±{m) and h±{m) operating on Rm{p) — Xm{p) / \fP-i 




(35) 
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When operating on the whole wave function, e*™''^i?^(p), a±(m) and 6±(m) may be re- 
placed by a± and b±, 

dp p 

a_ = e-'^[^ + -h-pl 
dp p 

h+ = e"t>[^--h-pl 
dp p 

6_ = e-"*'[^ + -h + p]. (36) 
dp p 

It is easily shown that, similar to (5) for a 2D hydrogen atom, for a 2D isotropic oscillator 
we have 

[Iz, a±] = ±a±, [Iz, b±] = ±6±, (37) 
i.e. both a± and b± are also angular momentum raising and lowering operators, but unlike 
the hydrogen atom, here a± and b± are not conserved quantities, 

[^,a±]=Ta±,[//,6±] = ±6±. (38) 

Therefore, the operators a± and b± themselves can not be directly equivalent to the 
conserved quantities characterizing the elliptic orbit. However, it is easily verified that 
the conserved quantities lz,Qi and Qxy can be constructed in terms of a± and b±, 

h = ^(a_a+ - 6+6_), 
Qi = -^(a_6_ + 6+a+), 

Qxy = -^{ci'-b- - b+a+). (39) 

V. SOME DISCUSSION ABOUT ID SYSTEMS 



The Schrodinger's factorization method and the concept of raising and lowering 
operators were extended extensively in supersymmetic quantum mechanics'^^"^^' to treat 
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the Schrodinger equation for a particle in a general ID potential V{x). It was shown that 
for a potential V{x), provided the ground bound state energy £'0 is finite (£'0 ^ —00) and 
the ground state wave function \l'o(.i') is differentiable, one can construct the corresponding 
raising and lowering operators, ^4+ and A, and the Schrodinger equation can always be 
factorized. The supersymmetric partner Hamiltonian, if_ = A^A and — AA^, have 
the same energy spectra, Ej^^ = -E^+i, (n = 0, 1, 2, • ■ • , ), except the ground state energy of 
{Eq~^ — 0), and the eigenstates with the same eigenvalue of H_ and are connected 
with each other by A'^ and A. It was shown that this is due to the shape invariance^"^'''! of 
V{x), which may also be considered as a special kind of dynamical symmetry 

It is interesting to note that the classical orbits of all bound particles in a regular 
ID potential are always closed. However, no classical close orbits exist for some singular 
potentials, e.g. , for a ID hydrogen atom (e = l)'^^] 

V{x) = — -j— r, { — 00 < X < +00). 

On the other hand, for a quantum ID hydrogen atom, the ground state energy Eq = —00, 



and "ifolx) ~ yS{x) is not differentiable at the origin. Thus, it seems understandable why 
the Schrodinger equation of a ID hydrogen atom can not be factorized. 

As for a ID harmonic oscillator, the energy eigenvalues and eigenstates are well 
known, which are quite similar to those for a 3D isotropic harmonic oscillator {I — 
case). The raising and lowering operators derived from factorization are = :^(^ ~ ^) 
and a = -^{x + ^), which connect the neighboring eigenstates with opposite parity 
(AA^ = 1). However, it should be noted that the ID harmonic oscillator potential formally 
corresponding to a 3D isotropic oscillator {V{r) = r^/2,r > 0) is 



whose energy levels are E^ — {N + 1/2), N — 1,3, 5, •. The usually adopted ID harmonic 
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oscillator is V{x) = (—00 < x < +00) with reflection symmetry, whose levels are 
En — N+1/2, N — 0,1, 2, 3, • and the neighboring eigenstates are of opposite parity. Prom 
this, one can understand why there exist two kinds of raising and lowering operators, A± 
and B±, for an nD (n > 2) isotropic harmonic oscillator, and A± and B± are different in 
form from the operators a"*" and a for a ID harmonic oscillator. However, one may use the 
product operator of A and B^^\ i.e. the operator C as the raising and lowering operators 
connecting the neighboring eigenstates with the same parity. In fact, for a 3D isotropic 
harmonic oscillator t^l 

C(, = 0,JVTT) = ^+r^-2rl-l 

is the same form as 2a+a — -£2 + x'^ — 2x^ — 1 for a ID harmonic oscillator, and both 
have the selection rule AN — 2. 

This work is supported by the National Natural Science Foundation of China and 
the Doctoral Program Foundation of Institute of Higher Education of China. 
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